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It is argued that the expected turn-down in a; — Q 2 of the cross sections (structure functions 
F2(x, Q 2 )), assumed to result from the saturation of parton densities in the nucleon, is re- 
lated to a phase transition from the (almost) ideal partonic gas, obeying Bjorken scaling, 
to a partonic "liquid". This can be quantified in the framework of statistical models, per- 
colation and other approaches to collective phenomena of the strongly interacting matter. 
Similarities and differences between the case of lepton-hadron, hadron-hadron and nuclear 
collisions are discussed. 



1 Introduction 

Based on different observations, models and equations governing deep inelastic scattering (DIS) 
and related processes, a "saturation" regime is expected when certain values of low enough x 
and relevant Q 2 are reached. According to the dipole model of DIS, this regime already has 
been achieved and it is characterized by the "saturation radius" [Tj i? 2 , = (x/xo) x /Qq, with 
Ql = 1 GcV 2 , x = 3 • 10" 4 and A = 0.29, found from a fit to the DIS data at x < 0.01. 
On more general grounds, saturation could be expected also from unitarity: the rapid (power- 
like) increase with 1/x of the structure functions/cross sections may suggest that unitarity 
corrections will tamper this rise, although formally the Froissart bound has never been proven 
for off-mass-shell particles, thus unitarity does not provide any rigorous limitation for such 
amplitudes [2 [3]. One more argument is physical: the rise of the structure function F2(x,Q 2 ) 
reflects the increase of the parton density (parton number in the nucleon). Since this number 
increases as a power, and the nucleon radius is known to increase as In s (or, at most In 2 s) 
(shrinkage of the cone), the particle number density within a nucleon increases, inevitably, 
reaching a critical value where the partons start to coalesce (overlap, recombine etc). The 
qualitative picture of this phenomenon in the Q 2 — 1/x plane is well known and cited in various 
contexts (see, e.g., [3]). 

Quantitatively, the dynamics depend on many, poorly known, details, such as the properties 
of the constituents and their interaction within the nucleon. 

Below we propose a novel approach to the saturation phenomenon in DIS and related pro- 
cesses based on the collective properties of the excited nucleon. Namely, we suggest that, below 
the saturation regeon, the nucleon in DIS is seen as a gas of almost free partons. With their 
increasing density, the constituent gradually overlap and, starting from a certain value of x and 
Q 2 , the gas of free partons coalesce condensing in a liquid of quarks and gluons. Saturation 
corresponds to the onset of the new phase. 
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Figure 1: Phase diagram of DIS, from Ref. 



The thermodynamic properties of such an excited nucleon are characterized by its tempera- 
ture, pressure etc, and a relevant equation of state (EoS). The statistical treatment of partonic 
distributions is by far not new, see e.g. [51(51 [7]. What is new in our approach, is the interpre- 
tation of the saturation in DIS as a manifestation of the transition from a dilute partonic gas to 
a liquid. The details (nature) of this (phase?) transition are not known. It can be of the first, 
of second order, or, moreover, be a smooth cross-over phenomenon. Our main argument is that 
the volume of the nucleon confining the partons (quarks and gluons) in the interior increases 
slower, at most as ~ In 6 s (more likely, as ~ In 3 s), while the volume occupied by the interior, 
quarks and gluons, increases as a power, thus resulting in a limiting behavior: a gas-to-liquid 
cross-over or a phase transition. The present contribution is a first step in understanding this 
complex processus. 

In a related paper, Ref. [5], the phase structure of the hadronic matter in terms of its tem- 
perature T and its baryochemical potential fi, was studied in the framework of the percolation 
theory. The percolation mechanism was used in [9] to obtain a limiting energy dependence of 
the hadronic matter at s — > oo . 



2 Saturation 

We define the saturation line (in the x — Q 2 plane) as the turning point (line) of the derivatives 



3(]nQ 2 ) ' xy ' (<91n(l/:r)) : 

called Bq or B x slopes, where F2(x, Q 2 ) is a "reasonable" model for the structure function, i.e. 
one satisfying the basic theoretical requirements, yet fitting the data. For example, the model 
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Figure 2: The surface of Bq(x,Q 2 ) calculated in Ref. jlOf . 



for F 2 (x,Q 2 ) of Ref. [10] interpolates between Regge behavior at small Q 2 and the solution of 
the DGLAP evolution equation at asymptotically large values of Q 2 , practically for all values 
of x. The resulting two-dimensional projection of the Q— slope is shown in Fig. 1. In our 
interpretation, the critical line (saturation=phase transition) occurs along the fold line on this 
figure (compare with a similar figure, Fig. 2 of Ref. [TTJ, derived from a different model). 

The main goal of this paper is the identification of this line(point) with the critical line(point) 
on the T, fi phase diagram of an excited nuclcon viewed as a thermodynamical system. The 
thcrmodynamical approach to DIS may provide a new insight to this complex phenomenon. 
We are aware of the limited time scales in a deep inelastic scattering from the point of view 
of thcrmalization, a familiar problem relevant to any thermodynamical description of hadronic 
systems. Let us only remind that the thermodynamic approach to high-energy scattering and 
multiple production, originated by Fermi and Landau's papers, were applied to hadrons, rather 
than heavy ions. 

3 Statistical models of parton distributions 

The statistical model of parton distributions in DIS was considered and developed in quite a 
number of papers [S[6j[7]. In its simplest version, one assumes [5] that inside the nucleon, the 
valence quarks, as well as the sea quarks and antiquarks and gluons form a noninteracting gas in 
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equilibrium. This simple picture may be further developed in two directions: 1. Introduction of 
the effect of the finite size (and its energy dependence!) of the nucleon on statistical expression 
for the number of states for the unit energy inerval; 2) Account for the Q 2 evolution, that can 
be calculated either from the DGLAP equation or phenomenologically, as e.g. in Ref.[TU]. A 
likely scenario emerging [SJE] for this high quark density n q = (n q — rig)/ 3 system is that quarks 
form Cooper pairs and new condensates develop. 

The nucleon of mass M consists of a gas of massless particles (quarks, antiquarks and 
gluons) in equilibrium at temperature T in a spherical volume V with radius R(s) increasing 
with squared c.m.s. energy s as Ins (or In 2 s). The invariant parton number density in phase 
space is given by [12] 

dn l dn = gf(E) 

d 3 p i d 3 rl d 3 pd 3 r ( 27r )3 ' I > 

where g is the degeneracy (g — 16 for gluons and g = 6 for q and q of a given flavor), E, p is the 
parton four-momentum and f(E) = ^cxp[/3(£' — fj,)] ± lj is the Fermi or Bosc distribution 

function with (3 = T~ x . Quantities in the infinite momentum frame (IMF) are labeled by the 
subscript i. 

The invariant parton density dn 1 /dx in the IMF is related to dn/ dE and f(E) in the proton 
rest frame as follows [5] 

dn* gV{s)M 2 x f M/2 



and the structure function 



dx (2tt) 2 

' dn 1 \ / dn 



dEf{E), (3) 

AI/2 



Without any account for the finite volume of the hadron, this SF disagrees with the data. Finite 
volume effects can be incorporated, following R.S. Bhalerao from Ref. [5], and the result is 

dn/dE = gf(E)(VE 2 /2tt 2 + aR 2 E + bR), (5) 

where V and R are energy dependent and a, b, in front of the surface and curvature terms, are 
unknown numerical coefficients. Their values are important for the final result, but they cannot 
be calculated from perturbative QCD. A rather general method to calculate these important 
parameters can be found in Ref. |13j . We intend to come back to this point in a subsequent 
publication. 



4 Percolation 

Percolation as a model of phase transition from colourless hadrons to a quark gluon plasma 
(QGP) was studied in a number of papers, recently in Ref. [8]. Below we apply the arguments 
of that paper to the saturation inside a nucleon, where extended (dressed) quark and gluons, 
rather then mesons and baryons considered in Ref. [5] percolate into a uniform new phase of 
matter. Both objects are coloured particle inside a colourless nucleon. We start with a short 
introduction to the subject. 

Consider TV spheres of radius i?o and hence volume Vb = (47r/3)i?g in a "box" of size V, 
with V > Vq. 
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Percolation (clustering) of spheres, in three spacial directions, was studied for the case of: 

a) arbitrary overlap [15] and b) for those with impenetrable hard core, allowing only partial 
overlap [5J. 

For the case a), the first percolation point (partial percolation), with 30% of occupation, 
occurs for the density n = N/V, n m sa 0.35/Vo, the largest cluster having the density of about 
1.2/Vq. A second percolation point occurs at n v = 1.22/Vb, when 70% of space is covered by 
spheres, i.e. the vacuum disappears as a large-scale entity. 

The existence of two percolation thresholds, one for the formation of the first spanning 
cluster of spheres and the the second one for the disappearing of a spanning vacuum "cluster", 
is a general feature of the 3-dimensional percolation theory. 

b) An impenetrable spherical core, with R c = Rq/2, and the spheres can only partially 
overlap. Here again one has two percolation thresholds, at h m m 0.34/Vb (close to the case a)), 
and vacuum percolation, at n v = 2.0/Vb, requiring a higher density compared to a). 

4.1 GPD H(x,t) and impact parameter parton distribution q(x,b) 

The number of constituents in a nucleon can be found by integrating in b (impact parameter) 
the general parton distribution, e.g. that of Ref. [14]. 



q{x,b) 




1 1 1 : > contrary to q(x), q(x,b) is dimensional, with a dimension of squared mass m 2 , inter- 
preted as the transverse size of the extended parton in the hadron, m = R~ x , q{x,b) = 
m 2 q(x 1 b), q(x,b) being the partons number density [2]. In paper [T3], m was a constant; 
here we choose it to depend on the the photon virtuality: m — ► mo ln(Q + m p ) = R^ 1 (Q), 
where m is the mass of the lightest vector meson, m = m p . Note the inequality q(x, b) < l/S q , 
where S q (Q) = ln _2 (Q + m p ). 

Thus, the nucleon is composed of N = 2ir f Q dx J °° bdbq(x, b, Q) extended partons with the 
transverse area S q (Q). 

4.2 A toy EoS 

One has for mesons 

Nm(T) = 3^T 3 . 

On the other hand, mesons percolate whenever n v = 1.22/V^, where Vh = (47r/3)i? 3 , and we 
use for the probe meson radius Rh = l.l/ln(Q + m p ), which in the limit of a real photon, 
Q — ► 0, matches the relevant value Rh = 0.8 fm used in Ref. [5]. Solving n h (T) = rif yields 

T M (T, Q) ~ 171 ln(Q + m p )MeV 

as the limiting temperature through meson fusion (cf. T n ~ 240 MeV of [5]). In a similar 
analysis of the phase diagram of hadronic matter [5], the (Q— independent) limiting temperature 
fx 240 MeV was obtained. 

N.B. The parameters appearing below should be rescaled with account for the replacement 
m — > ln(Q + m p )! 
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The density of point-like nuclcons of mass M is at T = 

7i b ( M ,T = 0) = ^(/i 2 -A/ 2 ) 3 
or, by using the van der Vaals approach of Ref. [5] 

n b (T,fi) 



n B (p, T = 0) 



l + n b (T, fi)V e 



With increasing nuclcon density, the empty vacuum disappears for n v ~ 2/Vh(Q) ~ 0.93 ln(Q + 
m p )) 3 /m -3 , which, for real photons (Q = 0) corresponds to about 5.5 times standard nu- 
clear density. Solving ub(T = 0,/i) = n v gives fi v ~ 1.12 ln(Q + m p )GeV for the limiting 
baryochemical potential at T = 0. 

In the region of low or intermediate \i in the Tfi diagram, one can approximate the density 
of point-like nucleons by the Boltzmann limit 



„ i(( ,,r)^(M)V 2 (M / r,e^^(f) I V-./, 



As a result, one obtains a family of curves for the phase diagram T(/z) and for the EoS p(T) 
for various values of Q and x (cf. Figs. 3 and 4 of Ref. [5]). 



5 Conclusions 

In this talk a new approach to the saturation phenomena in deeply virtual processes - DIS, 
DVCS, VMP - diffractive and non-diffractive - is suggested. The basic idea is a physical one: 
Bjorken scaling implies that the nucleoli in DIS and related processes is a system of weekly 
interacting partonic gas, that can be described by means of Bose-Einstein or Fermi statistics. 
As the density increases (with decreasing x and relevant Q 2 ), seen as the violation of Bjorken 
scaling, the system reaches a coalescence point where the gas condenses , eventually to a liquid. 
The thermodynamic properties of this transition can be only conjectured, and further studies 
are needed to quantify this phenomenon. 

Of two models presented in this note, more promising seems be the first one (Sec. 3). 
Further studies will show its viability. 

An alternative measure of the onset of saturation and expected change of phase can be 
related to non-linear evolution equations. Saturation and a phase transition are expected when 
the non-linear contribution overshoots the linear term. 

The phenomenon discussed in the present note may have much in common with the colour 
glass condensate proposed in the context of heavy ion collisions (see, e.g. [TB] and earlier 
references therein) . Apart from similarities (condensation of quarks and gluons as x — > 0) there 
are apparent differences: for example, hydrodynamical flow is not expected in DIS. In any case, 
the dynamics of the strong interaction is the same in lepton-hadron, hadron-hadron and heavy 
ion collisions. 
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